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The antiferromagnetic Heisenberg model on an anisotropic kagome lattice may be a good minimal 
model for real magnetic systems as well as a limit from which the isotropic case can be better 
understood. We therefore study the nearest-neighbor Heisenberg antiferromagnet on an anisotropic 
kagome lattice in a magnetic field. Such a system should be well described by weakly interacting 
spin chains, and we motivate a general form for the interaction by symmetry considerations and by 
perturbatively projecting out the inter-chain spins. In the spin 1/2 case, we find that the system 
exhibits a quantum phase transition from a ferrimagnetic ordered state to an XY ordered state as the 
field is increased. Finally, we discuss the appearance of magnetization plateaux in the ferrimagnetic 
phase. 



I. INTRODUCTION 

Frustrated antiferromagnets are of considerable inter- 
est because frustrating interactions lead to strong fluc- 
tuations both of classical (thermal) and quantum me- 
chanical origin. Unfortunately, even the simplest model 
Hamiltonians for such materials are often not easily 
analyzable by conventional theoretical techniques. An 
outstanding example is the nearest-neighbor quantum 
Heisenberg antiferromagnet on the kagome lattice. This 
structure is frustrated to a particularly high degree, and 
the extensive classical ground state degeneracy may be 
identified as a mechanism for strongly enhanced fluctua- 
tions. This strongly limits the usefulness of the standard 
semi-classical spin wave technique. The behavior of the 
spin 5* = 1/2 case in zero field is particularly puzzling, 
with different numerical ap proaches y ielding contradic- 
tory and/or puzzling resultJilSSHiniSEl Analytical stud- 
ies have not been any more illuminating as they necces- 
sarily involve numerous approximations that ofte n lower 
the symmetry of the problem substantiall y^'^'^"*^"'^ * ^^ * ^'^! 
However, it should be noted that recent experimental 
work on ZnCu3(OH)6Cl2 seems to fav or a gapless, mag- 
netically disordered ground stat d^^ * ^^ . 

One window into the behavior of frustrated antifer- 
romagnets is their magnetization process, i.e. the curve 
of magnetization versus applied field. Variation of the 
magnetization over the full range from zero to satura- 
tion opens another dimension of phase space in which 
to explore the phase diagram of these systems, and per- 
haps find physical and analytical insights. A qualitative 
feature to be understood is the occurrence of plateaux 
with quantized magnetization. The structure of these 
plateaux reveal some aspects of the correlations of the 
system. For the low-spin kagome lattice, there are sug- 
gestions of a plateau with M = 1 / SMg {Mg is the satura- 
tion magnetization) .ISESlIZIIIElIini This is a very natural 
structure for the kagome lattice, and can be understood 
as a state with two parallel and one anti-parallel spins per 
triangle. One experiment, however, suggests a plateau 
with AI = 1/2 A/3, which is not an obvious statd^D, 

In this work, we consider an antiferromagnet on the 
anisotropic kagome lattice illustrated in fig.[l] which may 




FIG. 1: The anisotropic kagome lattice with antiferromag- 
netic coupling J between nearest-neighbor chain spins {S 
spins) and coupling J' between S spins and inter-chain spins 
(/ spins). We make the approximation that the / spins are 
fully polarized by the external magnetic field. 



be viewed as a set of spin chains (whose spins we call 
S spins) coupled together through intermediate spins (/ 
spins). The anisotropy could stem from an inhomogene- 
ity in the spin interactions, or from an actual spatial 
anisotropy. A well known example of such a system is 
volborthite, in which spin 1 /2 moments residing on cop- 
per atoms form a kagome network that is anisotropic due 
to differing superexchange bond angles.f^^ Recently, there 
has been considerable progress in developing techniques 
to analyze such systems of quantum spin chains weakly 
coupled together by frustrating interactions. Such meth- 
ods have the advantage that they can apply directly to 
low spin and full Heisenberg (SU (2)) symmetry, and fully 
include the effects of quantum fluctuations. Here we de- 
termine the physics of the magnetization process in such 
a limit of the kagome lattice. 

We assume that the spins interact via a nearest neigh- 
bor Heisenberg interaction, and we take the magnetic 
field to be applied in the +z direction. Taking the con- 
vention that lowercase indices only sum over sites on the 
spin chains (S'-spin sites) while uppercase indices only 
sum over intermediate /-spin sites, we decompose the 
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Hamiltonian as 



II. PERTURBATIVE CALCULATION OF 
EFFECTIVE HAMILTONIAN 



H = 

Ho = 

H' = 



Ho + H' 



<»j> 



<Ai> 



Sill 



<Ai> 



Since we consider the limit J' / J <C 1, the low energy 
states are those in which the I spins are completely polar- 
ized by the external field, and so we work with an eS^ec- 
tive Hamiltonian projected into the space of such states. 
We approximately calculate this efi^ective Hamiltonian as 
follows. 



This decomposition is convenient in the limit we focus 
upon, namely J' ^ J and J' <C h. Note that we need 
not assume ft./ J is either large or small, and so can ex- 
plore (almost) the full magnetization process. The region 
around /i = is excluded however by J' <C h. We leave 
this more challenging limit for future work. We have cho- 
sen to incorporate the J' term coupling the z-components 
of chain spins to inter-chain spins into Hq, which is con- 
venient since that term commutes with the large field {h) 
term. This however means that the expansions in powers 
of H' and in J' do not strictly coincide, since the latter 
arise from two sources. 

Because J' is much smaller than the applied field, the 
system may be described by a low-energy effective the- 
ory involving only states in which the / spins are fully 
polarized. In the limit J'/ J 0, one expects that this 
low-energy effective theory will be a set of independent 
Luttinger liquids describing each chain (this will always 
be the case for 5 = 1/2 and will also hold for S = \ 
as long as the external field is larger than the zero-field 
gap to magnetic excitations). For very small J', we con- 
tinue to assume that the Luttinger liquid description is 
accurate, but there must also be weak inter-chain inter- 
actions. While the inter-chain interactions may be calcu- 
lated, in principle, by symmetry considerations alone, we 
find it convenient to use an approach combining pertur- 
bation theory and symmetry arguments that allows the 
couplings to be estimated by hand. 

For the case of a spin 1 system in a large enough field, 
we find that only one relevant interaction term emerges in 
an RG analysis of the general inter-chain coupling. Fur- 
ther analysis predicts that the resulting long-wavelength 
theory favors an XY antiferromagnetically ordered phase 
at zero temperature. 

For the case of a spin 1/2 system, however, an RG 
analysis of the inter-chain coupling reveals two relevant 
interaction terms. For large applied field, one interaction 
dominates the other in the long wavelength limit, leading 
to the same XY antiferromagnetically ordered phase as 
in the spin 1 case. However, for small applied field, the 
other relevant interaction dominates, leading to a fer- 
rimagnetically ordered phase. Finally, we reconcile the 
discrete translational symmetry of the microscopic lattice 
with our continuum theory and show that magnetization 
plateaux may appear in the ferrimagnetic phase close to 
the transition. 



A. Elimination of inter-chain spins 

Let P be the projection operator onto the subspace of 
all states with the / spins fully polarized, and let Q = 1 — 
P. Furthermore, for a general state let j-^i) = P\'(p) 
and IV'q) = Q\'ip) so that 1-0) = \^^) + IV'q)- 

Acting on the eigenvalue equation H\ij}) — E\il}) with 
P and then with Q gives the equations 



Letting R~{E — Hq)~^, ^ may be written as 
\^q) = RH'lij^) + RQH'li^Q) 

and then iterated to give 
1 



(1) 
(2) 



li'Q) 



1 - RQH 



-RH'ltp^) 



(3) 



Finally, upon replacing IV'q) in (IT|) by the expression 
in ([3]), we obtain an equation for alone: 

1 



Ho + PH' 



1 - RQH' 



RH' 



(4) 



Note, however, that Eq.Q is not really a linear eigen- 
value equation, since H^a itself depends on E through 
R. 



B. Second Order Approximation 

We would like to use the assumption that J' <^ J 
and J' <C ft to simplify the expression for H^tf and find 
an approximate eigenvalue equation for iV't)- First, we 
make the generalization 

H ^ Hc,= Ho + aH' 

in order to more conveniently keep track of powers of H' . 
Henceforth, take jV") to be the ground state of Ha. The 
ground state energy admits an expansion 

E — Eq + El + E2 -l- . . . 
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where En oc a" and Eq is the ground state energy of Hq. 
ff may then be expanded as 

with 



H2 = PH'RoH' , 



(5) 



and Ro = {Eo-Ho)-\ 

The second term in the expansion of H^a may be sim- 
plified by use of the fact that we are restricting it to act 
only on states in the image of P. One finds 




H2|Vt> 



FIG. 2: The generating symmetries of the space group cmm 
of the anisotropic kagome lattice: (a) translations, (b) ro- 
tations by TT about the hexagon centers, and (c) reflections 
through a vertical line passing through the hexagon centers. 



2 S 

<Aik> 



SI 



QZ 

nnA j 



iEo-Ho)-h + J'j: 
Within the P subspace, we may drop the ket and write 



H2 - y^H^, 



(6) 



with 



{aJ'fS 



is the Hamiltonian of a Luttingcr hquid on chain n. Here 
the "spin velocity" v and the Luttinger parameter g are 
known functions of h/ J and 5*. The variables 0„ and Lpn 
are dual boson fields living on the rfi^ chain and satisfying 

[LPn{x),9m{y)\ = -indnrnfix - y) . 

The physical interpretation of the phase fields is un- 
derstood from their relation to the microscopic spin op- 
erators: 



E SI 



<ik;A> 



iEo-Ho)-h+J'j: 



QZ 

j nnA j 



(7) 



where i, k are summed over the neighbors of A. Note 
that, because of the J' dependence of the denominator 
in Eq.(|7|, H2 contains non-trivial terms at all orders of 
J' greater than or equal to 0[(J')^]. This is convenient 
because it allows us to obtain some 0[( J')**] contributions 
by an only O(a^) calculation. 



III. CONTINUUM LIMIT AND SYMMETRY 
CONSIDERATIONS 

It is well-known that the low energy physics of an 
isolated antiferromagnetic Heisenberg chain in a uni- 
form applied field is described by Luttinger liquid theory, 
whenever the field range is such that the system remains 
gapless. In the limit J' ^ J, h, which we consider, all 
influences of the inter-chain spins upon the spin chains 
indeed occur at low energies, and so can be considered in 
this framework. 

Luttinger liquid theory consists of replacing the Hamil- 
tonian of each spin chain by that of a free boson field (and 
corrections that can be analyzed perturbatively) : 



Ho 



where 



-y 



1 



(8) 



(9) 



bi C0s(2ipn(Xj) ' 



]x,)] (10) 

Qx,) .(11) 



The coefficients ai, bo and bi depend on the chain mag- 
netization per site m^i, and, for S = 1/2, the wavevector 
Q is given by Q = 27r(moh -I- 1/2). 

From the above considerations, we expect that H2 (and 
higher order corrections) can also be expressed within the 
continuum (conformal) field theory. Formally, one may 
expand any local Hamiltonian in terms of the scaling op- 
erators of the decoupled fixed point theory. Specifically, 
the density can be written as 



nA,nA + 



lixA), 



(12) 



where O^^n+iix) is a local symmetry allowed operator 
involving degrees of freedom from the chains n,n + 1 
immediately above and below the site A, at horizontal 
position X, with scaling dimension A^-. This implies that 
the two point functions of these operators obey 



{o:,^n+iix)o:,,+,{x')) 



r'|2A„ 



(13) 



with constants that are dependent upon the conven- 
tion for normalizing the fields. The expectation value is 
taken in the continuum conformal field theory (CFT) de- 
scribing the decoupled Heisenberg chains in a field. The 
two point function of two difFerent operators ^ a' van- 
ishes if the operators have different symmetry (or more 
formally are descended from different primaries in the 
CFT). The sum over a may be thought of as being in 
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order of increasing scaling dimension. Terms with large 
scaling dimension are strongly irrelevant, and can there- 
fore be neglected. 

The first few terms in this expansion are strongly con- 
strained by symmetry. Considering the full microscopic 
Hamiltonian, the symmetries consist of the lattice space 
group {cmm - see fig.[2| and the U{1) rotational symme- 
try of the spins about the z axis. One finds 

n2 = Aj_ cos[6'„ - On+i] + K cos[2(^„ - 2i^„+i] (14) 

+ dxOn dx9n+l COs[6'„ — On+l] + . . . . 

Here we have suppressed the A subscript on n and x on 
the right hand side to keep the formula compact. 

The second order Hamiltonian, 7i^, is actually fur- 
ther constrained by an additional symmetry which one 
obtains only at this order. Specifically, for each term 
in associated with an inter-chain site A, the chains 
above and below this site can be independently reflected 
about a vertical axis through site A. This symmetry in 
fact requires that A_l = -I- 0[{J')^]. 



A. General prescription for coefficients 

In this subsection, we demonstrate that the above ef- 
fective continuum Hamiltonian is uniquely determined 
from the microscopic model, by showing how the coeffi- 
cients (A(j) may be obtained in principle from correlation 
functions of decoupled Heisenberg chains. In the follow- 
ing subsection, we will make some additional simplifica- 
tions in order to provide explicit expressions. 

In general, we can proceed by demanding equality be- 
tween Eq.(|7]) and Eq.([l4| (Eq.([l2|). The coeflicients can 



then be extracted by taking expectation values of these 
quantities with microscopic operators whose continuum 
operator content is known, i.e. 



(15) 



where ^ microscopic expression composed of 

lattice spin operators in the vicinity of site j (horizontal 
position x') on chains n,n + 1, and the omitted terms 
indicated by the ellipses contain only operators of larger 
scahng dimension than C5^.„+i. For our problem, the 

needed Q"^ operators and coeflicients are 



bl cos[6l„ -9n+l] 



oz nz 



1 

COs[2(/3„ - 2ipn+l] + -^dxfn dxifn+l + ■ 



— b^dxOn dxOn+1 cos[6'„ — 6'„+i] 



where ai,6o axe the coeflicients in Eqs. (10 111 
Knowing these, one has 



(16) 



iQ'j;'n.A,'nA + l'^2) 



Ca^<7{0n^n^_^_l{Xj)0'^^J^^^l{xA)) 
Cq-ActCct 



XA 



2A„ 



(17) 



The neglected terms decay faster with | I, and 

can therefore be distinguished from the dominant term 
above. Thus by calculating the left hand side in Eq.(17l, 
and extracting its long-distance behavior, one obtains the 
coefficient Ao-, since Ca,Ca are presumed known. 



From this formulation, we can already deduce the order of the various A^ coefficients. Consider flrst A^. We require: 



( ^tn A + i-nA^^^j,nA + l^^j + 1 , n a + 1 ) 



\x,^xa\'^'^ 



(18) 



The only obvious constraint to obtain a non-zero result for this correlation function is that the sites fc, i in Eq.j7| 
defining reside on neighboring chains ua , tla + 1 . There are 4 such combinations included in the sum in Eq. 
Taking the limit J' ^ in the resolvent denominator, we obtain 



{aJ'fS 



\Xn XA 



2A', 



i^jA,nA'^^jA + l,nA) 



1 



Eq — Hq — h 



(■^JA A + 1+ "^JA + 1 ,n A + 1 ) 



(19) 



There is no reason for the correlation function on the left hand side to vanish. To estimate it, we note that the 
resolvent denominator is negative for all eigenvalues and bounded below in magnitude by \Eq — Hq ~ h\ > h. Hence 
we estimate 



Al< 



(aJ') 



M2 



(20) 
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up to 0(1) coefRcients that can be smooth functions of h/ J . 

Consider now the remaining two non-zero coefficients, Az,A^. Both may be obtained from the expectation value, 
which we denote G: 



G 



qZ 



,{H^) = G2 + G3 + 0[(J')'] 



(21) 



Since the number of particles in chains ua and ua + 1 must be separately conserved, the sites i, A: in the sum in Eq.([7]) 
must be on the same chain to obtain a non-zero result. Hence we may write 



G 



{Ea-H^)-h+J'Y. 



QZ 



+ 

iA,nA 



S. 



(22) 



First consider G2, the second order in J' contribution obtained by taking the J' ^ limit in the above resolvent. If 
one takes this limit, then the only remaining operator on chain ua + 1 which is not translationally invariant is the 
single ^^^i- Applying a translation on chain ua + 1, then, one sees that the S'|„^_|_]^ can be replaced by the same 
operator at any other site of the chain, i.e. Sp „^4_i, with arbitrary j'. Hence, it can be replaced by average over 
all sites, i.e. the uniform magnetization, which is a good quantum number and non-fluctuation. Upon making this 
substitution, one finds 



G2 



< 



{aJ')^mch ( 5*. 



[SiA., 

™ch {Sj^uA [^j 



s 



1 



[^1 



-'jA + l.riA. 



+ 

jA,nA 



s. 



jA+l.riA 



(23) 



where the < is to be understood as indicating an inequality of magnitude (not sign) , using the same bound on the 
denominator as before. The final line is straightforwardly analyzed from known results for Heisenberg chains. Because 
there is only a single field (S'|„^) at position Xj, the long distance decay of G2 contains two pieces: an oscillatory 
term decaying with the exponent and a non-oscillatory term decaying with the exponent A^. The oscillatory term 
indicates a vanishing contribution in the continuum limit. The non-oscillating term has half the expected exponent, 
2A^ to correspond to a contribution to A'^. It instead represents the generation of a term proportional to dz^n, which 
generates a small smooth renormalization of the magnetization curve but is otherwise redundant. Thus neither A^ 
nor A^ are generated at 0[(J')^]. 

In order to calculate the leading contribution to Az and A^ then, one must calculate Gj, by expanding the resolvent 
in Eq.(22l to first order in J', giving 



G. 



< 




QZ QZ rc- 

^ j.UA j.riA + l I jA.r, 



QZ QZ rc- 

i,nA i,nA + l I j 



S. 



iA+l.MA J 



1 

(So -Ho)-h 



k nnA 



'^jtl + l,nA] 



S. 



jA+l,nAi 



k nnA 



s 



.jA+l,nA J 



2A^ 



TT^ \Xj 



- XA 



2A' 



(24) 



and so we find that Xz and A'^ are generally of order ( J') 



n3 



B. Explicit Calculation for Spin 1/2 System 



The above results provide a general formulation to 
calculate the coefficients in the effective Hamiltonian, 
Eq.( |14[ ), from correlation functions of a single Heisen- 
berg chain in a field. Because of the complicated form 
of 7i^, however, even these one-dimensional correlation 
functions are not simple to obtain analytically, or ex- 
tract from known results. In this subsection, therefore, 
we make some additional simplifying assumptions which 
allow an explicit calculation of the effective Hamiltonian 
for the spin 1 /2 system. This calculation confirms the 



general structure obtained above and allows for some 
semi-quantitative estimates. 

Consider again the the second order (in a) contribution 
to effective Hamiltonian density 



aJ' 



<ik-.A> 



iEo-Ho)-h+J'j: 



QZ 

j nnA '-'j 



s 



When the operator above acts within the low energy sub- 
space, one would like to say that both the term propor- 
tional to J' and the {Eq — Hq) term in the denominator 
are small, but it is not clear that {Eq — Hq) is small 
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unless it acts directly on the low energy states. So, we 
exchange the operator and the resolvent using spin 
1 /2 anticommutation relations to obtain 



aJ' 



where 



Rai-.o 



<ik;A> 



{Eq — Ho) 



j nnA 



-2hS^ + 2jJ2iS-S] + S:[SJ) 

j nni 

However, to be able to expand the RAi-a operators, 
we must make two rather unphysical assumptions. First, 
we take J <^ h and second, we only let the magnetic 
field couple to / spins. This second assumption may not 
be too bad, however, since our eventual replacement of 
the chains of S spins with Luttinger liquids takes the 
effect of the external field into account via the Luttinger 
parameter. 

After making the approximations just discussed and 
expanding, we finally obtain the following leading contri- 
butions: 



{aJ'f 



Ah 



<ik;A> 
l\2 11 



E ^t^] 



<ijk;A> 



{aJ'fj 

4/i2 



E E "^fe ^t^j Sf, s^s^ 



<ik:A> j iin i 



(25) 



In what follows, we set a back to 1. 

To connect this result to the continuum limit of the 



interchain coupling Eq.(14l, we now proceed to bosonize 



the spin operators through the identifications in Eq.(lO) 
and Eq.(lll. Keeping only contributions that involve 



products of operators from both chains, we find from the 
first term 



E s^s^ 

<ik-.A> 



JA-riA 
+ [UA UA+l] 



l) 



which confirms that is of order J'^, as we already 
argued by symmetry. 



The second term in Eq.(25 1 must be treated with more 



care, however. Now, we only need to keep operators that 
conserve particle number on each chain, since we have 



already found the leading contribution to and since 
A_L = 0[(J')^] by symmetry. After also dropping opera- 
tors that do not have contributions from both chains, we 
find 

<iki\A> 

^iSjA,nA + l + SjA + l,nA + l) + i^A ^UA+l] 
= (^jA,nA + l 



+ {SiA. 



Q+ 

riA jA + i^nA 



-^t ,nA SjA + l,nA})+ ["^ ^ "^^+1] ■ 



It is the last operator in curly braces (the dimerization 
operator) that must be bosonized carefully. Since it in- 
volves products of operators at very short distances, all 



subleading bosonic operators omitted in Eq.( 10 1 must be 
resummed, leading to the following continuum limii!22l 



CQ+ — dxipn-Cl C0s[2(pn+Qx]. 
IT 

(26) 

After using the above identity and collecting all relevant 
terms, we find that 

^ 5-5+5| = ^(c^-4)a,^„,a,^„,+i (27) 

<iki;A> 

- {2a\ + aici)(l + cos Q) cos[2(^„^ - 2(^„^+i] + . . . 

demonstrating that \z and A'^ are both of order J'^ . The 
last term in Eq. ( 25 1 does not generate any relevant in- 



teractions in the bosonic theory. 

To summarize, then, we have explicitly calculated the 
couplings in the effective Hamiltonian Eq.(14l for the 



spin 1/2 system up to third order in J' and the re- 
sults confirm the predictions of our symmetry arguments 
(which apply to systems with any spin). The explicit 
couplings are 

A± - 
Al = 

A. - 

At = 





2h 

2h 



,{2al 



aici)(l -l-cosQ) 
(4-4) . 



(28) 



IV. ANALYSIS OF THE CONTINUUM MODEL 

We have arrived, then, at a low energy effective the- 
ory describing the Heisenberg antiferromagnet on the 
anisotropic kagome lattice that takes the general form 



E 



Ztt q Ztt 



A_L C0S[6'„ - On+l] + Xz COs[2<y9„ - 2lfin+l] 



Aj^ dxOn dxOn+1 CQs[6n — On+l] 



(29) 
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A. RG Analysis 



the solutions to the flow equations for Aj^ and Xz are 



To better understand the behavior of the above model 
at large length scales, we employ the momentum space 
RG procedure, lowering the cutoff from A ^ A — dK 
and defining dl = dA/A. Including the marginal term 
X'z dx^Pn dxfn+i in the fixed-point action and expanding 
all quantities to first order in leads to the following 
tree-level flow equations for the other couplings (note 
that X'z may be regarded as a constant independent of 
the scale of the effective theory) : 



dX'^ 
~df 
dXj_ 
~df 
dXz 

~dr 



1 

2-2.g 



ttA' 

+ ^ 



ttA^A^ 
16u 



A', 



(30) 



These flow equations imply that for A^ 7^ and 
A^ 7^ 0, a finite value of Aj^ will be generated in the 
long-wavelength theory. Also, one can see that the 9- 
interaction (with coupling Aj_) is relevant for g>\ — 
while the (/^-interaction (with coupling A^) is relevant for 

g < 1 — 

Now for the case of spin 1 chain, Fath has determined 
numerically that g is always greater than one (see FatlPl, 
fig. 5). In particular, g is one at a chain magnetization 
per site of m^i, = 0, reaches a maximum value of 5 « 1.46 
at rn^h w 0.36 and approaches one again in the limit of 
full polarization. So, for a spin 1 system described by 
(29 1, the 6'-interaction is always relevant, while the ip- 



interaction is always irrelevant. 

On the other hand, for a spin 1/2 chain, the exact 
Bethe ansatz solution shows that g increases nearly lin- 
early from 1/2 to 1 as the chain magnetization per site 
increases from m^i, = to m^i, = 1/2 (see Hikihara and 
FurusakP^ fig. 2). So for a spin 1/2 system described by 
( 29 1 , the (/^-interaction and the ^-interaction are nearly 
always relevant, and we must find a further criterion to 
determine which interaction dominates the behavior of 
the system. So, we solve the fiow equations and deter- 
mine which coupling constant becomes ^ 1 first as one 
continues to integrate out high momentum modes. 



B. Competing Interactions in the Spin 1/2 System 



We take the initial conditions given by Eq.(28 1. Then, 
defining 



/3± = (^2- 
/3. = ("2-23- 



1 ttA 

2g ^ ^ y ' 

1 -kX'. 

1 ' 

2g Av 



Aj.(0 



J'5 A^bl 



0±l 



/i3 647r 

J'3 



^(4-c;,) C 

V V 
(4-c;,) e^-' 



XS) = ^(aiCi-f-2a2)(l + cosQ) e'^^' . 

We can ignore the term '--^ e^^^ in A_l(Z) since it rapidly 
goes to zero unless g>l/X'z- Identifying e' as L/a, i.e. 
the ratio of the renormalized cutoff length scale to the 
bare lattice scale, we find that 



A±(/) 
A.(/) 



J'- 



If we then define L±^ and as the length scales at 
which the 9 and Lp couplings, respectively, become ^ 1, 
we find that 



{J') 



/^5//3_L-3//3, 



The exponent is positive for Q < g < gc and negative for 

9c < 3 < 1, where gc ~ 0.636 — 0.635^ -I- . . . This means 
that for g < g^, ^ L±_ and the (/? interaction is the 
dominant one. For g > gc, the 9 interaction is dominant. 



C. Ordered Phases in the Gaussian Approximation 



1. Ferrimagnetic Phase 



First, consider a system described by the Hamiltonian 



(29 1 in the limit where g < gc, that is a spin 1/2 system 
with magnetization per site between to = and to w 1/5 
(on a scale where full polarization is to — 1/2). In this 
limit, the dominant relevant coupling is the (/j- interaction 
with coupling Xz , so we drop all other couplings not in the 
fixed-point action. The 9 fields may then be integrated 
out to yield an action entirely in terms of the ip fields: 



S 



E 



V 
2TUj 



+ X'zdx^Pndx^n+l + Xz C0s[2(^„ - 2(p„+i 



(31) 



Since the interaction term in the above theory is rele- 
vant, the long- wavelength modes of the ip fields will fluc- 
tuate very little from chain to chain if the coupling Xz is 
negative, while they will differ by 7r/2 from chain to chain 
if the coupling is positive. If Xz is positive, then, redefine 
the (/? fields by ipn{x) ipn{x) + ^n. The coupling term 
for these new fields will now be negative. 

After the appropriate redefinition, fluctuations of the 
ip flelds from chain to chain will be suppressed at zero 
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FIG. 3: An example of a ferrimagnetically ordered ground 
state of the spin 1/2 system with magnetization m< 1/5 and 
< 0. 




FIG. 4: The XY antiferromagnetically ordered ground state 
that occurs for both the spin 1/2 and spin 1 systems in a large 
enough field (Ax < case shown). 



temperature, and so it is reasonable to expand the cosine 
term in the action to second order about zero. Upon 
performing a fourier transform in the chain index n, the 
resulting gaussian action is 



2TTg 



dk 



X,T 



(1 



A3Cosfc)|9^(/5fe| 



1 



\drVk\ 



that S^{x) is a continuum operator with all fluctuations 
on length scales below integrated out. Since g < 1 in 
the spin 1/2 model, this factor will tend to suppress the 
oscillations in the magnetization. 



(32) 2. XY Anttferromagnetic Phase 



where the effective mass is defined as 
^(l-cosA:)A^ . 



(33) 



A system described by the action ( 32 1 has excitations 
with a dispersion relation 

^lip,k) = (1 - ^X'^cosk)p\^ + ml{k)v' (34) 

which vanishes as p, k — > 0, implying the existence of 
gapless excitations. Indeed, the action (32) is invariant 
under the transformation 

ipk=oix) ^ ^Pk=o{x) + (p {ip constant) 

because of the invariance of the system under translations 
along the chain or x-direction. The gapless excitations 
are the Goldstone modes associated with the breaking of 
this symmetry. 

Restoring the original definition of the fields </3„(x) as 
appropriate, we see that for small enough g, the ipn fields 
will be ordered as {ipn{x)) = + f n if Az is positive and 
{(fn{x)) = (p a Xz is negative. Then, from Eq.( 11 1 we see 



that the spins will be ordered ferrimagnetically: 
TO,h + (5rn(-l)"cos(2(^ + (3a;) if > 



TO^h + 6m cos{2ip + Qx) 



if Az < 



(35) 

(see fig. [3] for an illustration). The amplitude of 

the ferrimagnetic spin density wave oscillation, 5m ~ 



Now, consider a system described by the Hamiltonian 



(29 1 but in the limit where g > gc, namely a spin 1/2 



system with magnetization per site m > 1/5 or a spin 1 
system in a sufficiently large field (larger than both the 
Haldane gap and J'). In this limit, the dominant relevant 
coupling is the ^-interaction with coupling Aj_ , so we drop 
all other couplings not in the fixed-point action. 

We proceed much as before, integrating out the ip fields 
and redefining the 9 fields if Aj^ is positive. In this 
case, however, the appropriate redefinition is 9n{x) — > 
9n{x) ~\- Tin. Once the fields are defined such that the co- 
sine interaction makes them slowly varying at zero tem- 
perature, we expand the cosine to obtain a gaussian the- 
ory once again, but with an effective mass 



ml{k) 



2tt 



(1 — COS k) Aj 



(36) 



and a dispersion 



ujI{p^ fc) = (1 - ^A; cos k) [p2„2 + mlik) w^] . (37) 

The gapless modes in this limit are the Goldstone 
modes associated with the breaking of the shift symmetry 

9k=o{x) — > 9k=o{x) + 9 {9 constant) 

associated with the global rotational symmetry of all 
spins about the magnetic field direction. 

Restoring the original fields, we expect that, for large 
enough g, {9n{x)) = 9 + nn if \± > or (0„(x)) = ^ if 
Aj^ < 0. From the identification ( 10 1, we therefore expect 
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an XY antiferromagnetic ordering of the spins: 



j+n 



(-1) 



if A_L > 
if A_L < 



(38) 



where mg ^ bQ{J' / J)^^''^^^^^ is the staggered magneti- 
zation, which is reduced by fluctuations below the scale 
Lj_. See fig. |4] for an illustration. 

D. Lattice Effects in the Spin 1/2 Ferrimagnetic 
Phase 

1. Commensurate Ferrimagnetic Order 

For a spin 1/2 system with g < gc, the ordering of 
the zero mode of the ip field corresponds to the break- 
ing of the translational symmetry of the system in the 
X direction (along the chains). However, this symme- 
try is not actually continuous in the microscopic theory, 
but reduces to the discrete translational symmetry of the 
lattice. 

In order to reconcile the continuum and microscopic 



theories then, we observe that the Hamiltonian (29 1 may 



also contain terms of the form cos[2p (y9„(x)] for certain 
values of p that must be determined by symmetry. 

Under translation by one lattice spacing along the 
chain direction, the fields transform as 



(fin{x) ^n{x - 1) - (3/2 



(39) 



So, for a term of the form cos[2p (/3„(a;)] to respect the 
translational symmetry of the system, it must be that 
pQ — 27rfc, where k is an integer. Moreover, the action 
should be invariant under ipn{x) —> ipn{x) -I- tt, so fc must 
be chosen to make p an integer. 

Recalling that Q — 27r(m<,i, -I- 1/2), both requirements 
imply that only the set {pj} are allowed, where 



(40) 



and kj is the jth positive integer (in increasing order) 
such that Pj is an integer. We may identify 



A = 



27r 



(41) 



as the wavelength of the (continuum) ferrimagnetic wave 
(35 1. But, for a general magnetization, the wavelength 



A will not be an integer multiple of the lattice spacing. 
Instead, the ferrimagnetic wave above, when restricted to 
values of x that coincide with the lattice, will appear to 
be periodic with an effective wavelength A„ff = fco-^ = Po- 
So then, given that for a fixed magnetization the cor- 
responding term cos[2po ^n] will appear in the Hamilto- 
nian, it will bring the discrete translational symmetry of 
the lattice to bear on the ordered phase by restricting the 
values of {(pn{x)) = if to coincide with its minima. For 




Ferri- 
Magnet 



XY 

Antiferromagnet 



Saturated 



t 



FIG. 5: Schematic graph of the magnetization curve and 
phases of the spin-1/2 system. Note the hypothetical plateau 
just below the crossover point. The shaded region indicates 
the parameter space where our model does not apply. 



a negative coefficient of cos [2po this will force the 
ferrimagnetic oscillations to be symmetric about a par- 
ticular lattice site. On the other hand, if the coefficient 
is positive the oscillations will be symmetric about the 
point on the lattice between two particular neighboring 
sites. 



2. Magnetization Plateaux 

In addition to pinning the ferrimagnetic order to the 
lattice, the symmetry-allowed cos[2p (^„(a;)] terms in the 
Hamiltonian may lead to the appearance of plateaux in 
the magnetization curve of the spin 1/2 system. 

For general (but small) values of the external field, the 
ground state of the spin 1/2 system will be, in bosonic 
terms, a state in which the (pn{x) fields on each chain fluc- 
tuate about some constant average value as we discussed 
in detail above (we take the coupling Xz to be negative in 
this section for simplicity - the other case follows straight- 
forwardly). In such states, the chain magnetization per 
site smoothly oscillates as a function of x about its av- 
erage (to zeroth order in J') value m^,,. Note that to^,, 
is simply the magnetization of a decoupled Heisenberg 
chain and as such, is a monotonially increasing function 
of the external magnetic field 

Now consider some particular magnetization 
(which corresponds to a chain magnetization to£ — 
3m^/2 — 1/4) and first assume that the external field 
has been tuned to a value such that m^i-,{h^) = m£. 
In such an external field, the most important symmetry 
allowed term of the type cos[2piy9„(a::)] that may be added 
to the Hamiltonian is the one with the smallest value of 
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p given by Eq. ( 40 ) , namely 
fco 



Po 



(42) 



where ko is defined as the smallest positive integer such 
that Po is an integer. 

To illustrate, then, how a plateau may arise at this 
magnetization, we consider the Hamiltonian 



1 



+ 1 



- |A^o|cos[2((^„ - (Pn+l)] - ?7oCOs[2po'/'n + 6x] + ... 

(43) 

and no longer think of the external field as fixed, but 
as varying in some small range about . Since the co- 
efficient 5 is zero for h = by the definition of p^, 
symmetry considerations show that it must be given in 
general by 



S{h) = 2Trpo[m,i,{h) - m£] 



(44) 



as the momentum cutoff is lowered to l/L^ under the 
RG procedure. Since X^o is small (generally ~ J'^) and 
g < I for the spin 1/2 system in its ferrimagnetic phase, 
this renormalized value of 77 will decrease rapidly as a 
function of pq . 

Finally, then, we determine if a plateau develops 
around the magnetization by considering the com- 
petition between the kinetic energy terms in H^ft and the 
term with coupling rj. The latter term is minimized if 
{ipn{x)) = ~2jH^ ^ (that is, if {S^{x)) = m£ regardless 
of the external field - a plateau). However, such a field 
configuration costs some kinetic energy and the resulting 
total change in the linear energy density is 



(48) 



Therefore, we see that this arrangment remains favorable 
(i.e. the plateau exists) as long as 



Now, if we integrate out the short wavelength modes 
of the system in an iterative RG procedure (but with- 
out rescaling the high momentum cutoff), the effective 
coupling will vary as 



A,[i] - A,o (i) 



-2g 



(45) 



until we reach a length scale defined such that 
Az[iz] ^ ^/{Lz)"^ (i.e. this is the scale such that A^ would 
be 0(1) if we were rescaling at every RG step). There- 
fore, we find that 



(A.o) 



-l/(2-2g) 



(46) 



At this scale, the chains will be strongly coupled and the 
(p fields will vary slowly from chain to chain. 

Then, since the ip fields vary slowly from chain to 
chain, it is reasonable to take a continuum limit (pnix) — > 
ip{x,y)/Vd (where d is the chain spacing). Within 
such an effective (2 -|- l)-dimcnsional theory, the long- 
wavelength ffuctuations of the field (p{x,y) are bounded 
and so the term in the action rjo cos[2pQ ipn], which may 
have been an irrelevant perturbation to a decoupled chain 
is now a relevant perturbation, regardless of the value of 
Po- 

Now, since the scaling dimension of the cos[2po</3] term 
is Pq g, the effective coupling 77 becomes 



77(L,) ^ ,70 • (L,)-foff 



'70 • (Azo)"-"« 



(47) 



\5{h)\ < ^2pQ 



2Tig 



(49) 



|m,„(/l)-TO£| < (j')3Po9/4(l-ff) _ 



A significant plateau will only appear, then, for a rel- 
atively small which is itself determined by m^^ 
through (|42|. 



We therefore expect the widest plateau, if one appears 
at all, to appear just below the crossover point from the 
ferrimagnetic to the XY antifcrromagnetic phase in the 
spin 1/2 system (see fig. [5]). Furthermore, we may esti- 
mate the dependence of the plateau width on J' from our 
earlier estimate that gc ~ 0.636 and from the fact that 



the smallest po allowed by Eq.(42l near the crossover is 
about 9. Using these values in Eq.(50l, the power of J' is 
^ 100, so any plateau near the crossover will be narrow 
indeed. 
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